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Speeding protein folding beyond the Go model: How a little 

frustration sometimes helps. 

Steven S. Plotkin 
Department of Physics, University of California, San Diego 



Abstract 



Perturbing a Go model towards a realistic protein 
Hamiltonian by adding non-native interactions, we find 
that the folding rate is in general enhanced as ruggedness 
is initially increased, as long as the protein is sufficiently 
large and flexible. Eventually the rate drops rapidly to- 
wards zero when ruggedness significantly slows conforma- 
tional transitions. Energy landscape arguments for ther- 
modynamics and kinetics are coupled with a treatment of 
non-native collapse to elucidate this effect. 



I. INTRODUCTION 

Theorists seek to capture the essence of protein 
folding with simple models of a self-interacting poly- 
mer chain §|,|,HH,||,!|,|!,!|0. Th<3re are 
two distinct limits pertaining to the nature of the in- 
teractions in this minimalist approach. One is that 
of purely random interactions, and is considered too 
frustrated to describe real proteins. Another is the 
Go model p7|, where the polymer is self-attractive 
only for those parts of it in their native configura- 
tions. This is considered too unfrustrated to describe 
real proteins, and also impossible to achieve in prac- 
tice. As these two models bracket the behavior of real 
proteins, we consider perturbing from the Go model 
towards real protein interactions by adding some non- 
native heterogeneity. Some of the effects of adding 
non-native interactions on the folding mechanism for 
the Honeycutt-Thirumalai /3-barrel model ]17[ were 
investigated in |2^,g2[. At first glance one would ex- 
pect that adding frustration begins to slow the rate at 
the transition temperature, or at best has initially no 
effect. What follows is a derivation of the somewhat 
counterintuitive result that in general the folding rate 
initially increases as ruggedness is increased from zero. 
Eventually of course the rate decreases drastically, so 
a plot of the folding rate vs. the amount of non-native 
heterogeneity should look like fig. ([!]) . Then the ques- 
tion of where real protein interactions reside on this 
plot may be addressed. For some fast-folding proteins, 
it is possible that non-native noise in the system may 
actually assist folding. 



II. THERMODYNAMICS 

Consider first the thermodynamics of a protein ob- 
tained from a statistical analysis of a correlated land- 
scape p9[ |. The energy, entropy, and free energy as 
functions of the fraction of native contacts Q, are 
given by [] 



E(Q) = QE N - 
S(Q) = S C (Q) - 



A 2 (Q) 
T 
A 2 (Q) 
2T 2 



(1-0) 



(1-0) 



F(Q) = QE N - TS C (Q) - 



A 2 (Q) 
2T 



(i-O) 



(la) 
(lb) 
(1c) 



These quantities are shown in figure (0), and the pa- 
rameters used in them are given in table 1. S C (Q) in 
eq. dig) is the configurational entropy in the system 
vs. Q, E N is the extra internal energy in the na- 
tive state (the stability gap), and A 2 (Q)(1 — Q) is the 
non-native variance which is a measure of the overall 
ruggedness of the energy landscape (see eq. ||). The 
variance A 2 (Q)(1 — Q) — ► as Q — ► 1. In obtaining 
the functional form of the non-native ruggedness, it is 
assumed here that all the native contacts have roughly 
the same strength, g 

The native energy is the number of native contacts 
M times the mean native attraction energy e (e < 
0). If N is the number of interacting residues in the 
polymer chain and z is the number of effective bonds 
per residue, 



E N = Me = zNe . 



(2) 



The scale for the overall non-native ruggedness 
A 2 (Q)(1-Q) is given by 



A 2 (Q) = Mb 2 V (Q) 



(3) 



1 We will generally set Boltzmann's constant k B = 1 in 
this paper, so temperatures have units of energy, and en- 
tropies are in units of k B . 

2 When there is variance in the native energies, the non- 
native ruggedness terms are proportional to (A 2 (Q) + 
QA„(Q))(1 — Q), where Aj is the native variance ]30| ]. 
If the set of native energies has variance but the distribu- 
tion is fully specified, then the ruggedness terms are again 
proportional to A 2 (Q)(1 — Q) feq|. 



where r)(Q) is the non-native packing density (0 < 
i](Q) < 1), b 2 is the intrinsic variance per non-native 
interaction, and M is the total possible number of 
(non-native) interactions, i.e. the native state is as- 
sumed to be fully collapsed with the maximal num- 
ber of contacts, and this is the maximal number of 
total contacts of any state. The de nsity r](Q) tends 
to increase upon folding (see section [I A), hence the 



ruggedness scale |A| increases as well. The strength b 
of non-native interactions is taken to be weak: 



where Arft- = rj(Q^ = 1/2) — 7/(0) is the change in 
non-native density between the barrier peak and un- 
folded state (Ar^ < 1). So long as Ar]^ > 0, the bar- 
rier height decreases as non- native heterogeneity (b 2 ) 
increases, as shown in fig. (g). We now show that this 
is nearly always the case, by considering the physics 
of collapse for our problem in question. 



A. The Collapse Transition 



e 



(4) 



therefore the ratio of folding transition temperature 
T F to thermodynamic glass temperature T G is large 






» 1 



(5) 



i.e. the proteins we consider are strongly (but not 
infinitely) unfrustrated- we are perturbing away from 
the Go model. 

The configurational entropy S C (Q) has the property 
that entropy loss on folding is more rapid initially than 
in later stages. We approximate this effect here by 
assuming the form 



S C {Q) = S (l - Q) - Tent(Q) 



(6) 



where S = Ns is the total conformational entropy 
in the unfolded (Q = 0) state (s is the log number of 
conformational states per residue), and Tent(Q) is a 
tent function: 



In this section we investigate the coupling of non- 
native density with the amount of native structure 
present in a protein, by showing that native topologi- 
cal constraints can induce a collapse transition on the 
non-native parts of the protein. Then the trend in 
eq. (M) of adding non-native ruggedness would be to 
lower the folding barrier. 

Collapse occurs below a temperature Tg, defined as 
the temperature where the free energy of the coil and 
collapsed molten globule phases (both at Q = 0) are 
equal: 



F co a{Te) = F mg (Tg) 



(11) 



Again using eq. (ftq), but now noting that the confor- 
mational entropies are different in the coil and globule 
phases, and that rj = 1 in the globule and r\ = in 
the coil phase, we have 



TqSc 



~TgS„ 



Mb 2 
~2Te~ 



Ma. 



(12) 



Tent(Q) 



2(f>Q 



Q<l/2 



20(1 -Q) Q>l/2 



Note we have now allowed for a mean homopolymer 
(,'J attraction a in general, for reasons which will become 
clear below. Using the reduction in entropy for col- 
We've let the barrier be at Q* = 1/2 for simplicity of la P sed vs - coil chain statistics |,|| 
argument. 

At the transition temperature T F , the free energy of Scoil ~ Sm 9 = N lo § v ~ N lo &\ v l e ) = N (I 3 ) 

the folded and unfolded states are equal: 



gives for the collapse temperature 



F(0)^F(1) 

A 2 (0) „ 



(8) 



Using this relation in eq. ([k]) gives 
F(Q) - F(0) 




(14) 



T 



(9) 



T F 



Note from (Il4) that when b = 0, Tg = za, i.e. the col- 
lapse temperature is the mean homopolymer attrac- 

— T t(n\ (1 — Q) r /^j\ /(-.-tion times the number of contacts per residue, and 

- \vn\UJ) 2T z [Tj\H) - V(V^ en a = QTg = by/7/2, i.e. non-native hetero- 

geneity can drive collapse, with the collapse tempera- 
ture now scaling with the root number of contacts per 
residue times the width of interactions. 

Now we note that in our model (Go perturbed by 
weak heterogeneity, with a = 0) collapse and folding 
will tend to occur together, with folding driving the 
collapse through native structural constraints. So the 
total density increases from zero to one as Q — > 1, and 
Arj , (10) the non-native density Tj(Q) should increase as well 



When 6 = the free energy at T F is the tent function 
(see fig. |), and so cj> in eq. (0) is thus F^(b = 0)/T F . 
Then from eq. (M) the free energy barrier at T F is given 
by 

AF* AF^(b = 0) Mb 2 



T F 



4T< 



since non-native polymer is more strongly constrained 
by larger native cores, see figure pi The simplest ap- 
proximation to capture this increase in density upon 
folding is to replace the mean homopolymer field a by 
the native energy scale e times the fraction of native 
bonds made Q: 



a{Q) = eQ . 



(15) 



This is the effective homopolymer field for the ensem- 
ble of states with fraction Q of native structure. Us- 
ing (fl3) in (nj) and noting that the glass temperature 



T a 



zb 2 



2s 



(16) 



my 



and stiffness effects. Applying bulk thermodynamics 
to the residual segments of non-native polymer may 
not be an accurate approximation in some cases. 

More complete treatments of the coupling of density 
with native similarity can be made within the energy 
landscape framework 29 1. It is fairly straightforward 
to write an approximate free energy as a function of 
both i] and Q and then minimize with respect to r\ to 
obtain the density as a function of Q. We have taken 
the simplest approach here to illustrate the coupling 
of collapse with the thermodynamics. Some caution- 
ary notes are made in Appendix A regarding a pos- 
sible reversal of the trend on barrier height in small, 
stiff proteins, or proteins with a significant amount of 
generic self-attraction. 



gives 




^=T F Rr + i = +«».[£■ ("0 






where s and s mg are the entropy per residue in 
the coil and globule state respectively. Note that in 
eq. (Tl7|) Tg > T F as long as the term in parenthe- 
ses is greater than one. This gives a critical value 
Qg where collapse occurs during folding, i.e. when 
Q ~ Qe, V ~ 1 and when Q £ Qg, r\ « 0. This is 
sketched in fig (f|A) below. Solving for Qg gives 



Qe 



1 



So V T F 



T, 



1 



s n b 2 



1 



s 2ze 2 



Tg(Q) = Qs a T F = Qze 



(18) 
(19) 



since, by construction of the problem, eqn. (j4|) holds, 
e.g. say b/e is about 1/20. Then the second term 
in (|l8|) is of order 1/400 and can be neglected. Equa- 
tion Q18D says that the more chain entropy the poly- 
mer has (the more flexible it is) the sooner it collapses 
when folding at the transition temperature. 

Calorimetric measurements of the conformational 
entropy change per residue in unfolding to the coil 
state for say barnase give s = 55 J/K • mol residue = 
6.8fc B per residue |20|| . This entropy also counts side 
chain conformational entropy, which is estimated to be 
about 13 J/K • mol residue = 1.6 k B per residue pOfl , 
giving a net chain conformational entropy of about 
5.2 k B per residue in the coil state, and therefore 
Qg = 0.2. For typical off-lattice simulations [p|, p3| 
s = 3.4fc B , therefore Qg = 0.3. So collapse typically 
occurs before the barrier is reached (see fig. (||)C). 
For lattice simulations Qg = 0.6, which is around 
Q^. In any event, T)(Q^) will tend to be greater than 
77(0) as long as the system is large enough and bulk 
thermodynamics can be used (however see caveats in 
appendix M), see figure 1. The values of non-native 
packing density obtained from simulations appear to 
be smaller than one, probably because of finite size 



III. KINETICS 

What about the rate? The question is now whether 
the increase in prefactor is larger than the decrease 
in barrier height, as we add non-native heterogeneity. 
Since the ruggedness is weak, the kinetics are single 
exponential (there is a single dominant folding bar- 
rier), and a Kramers law holds for the rate: 



k 



- l (b)e 



-AF*{b)/T 



(20) 



where the prefactor is proportional to the reciprocal 
of the reconfiguration time scale |3|, po[ |35|, |39J1 . 

If we were to follow the argument for the depen- 
dency of the prefactor on ruggedness for an uncor- 
related landscape pi, or for a correlated landscape 
at low temperature with activated dynamics |B9|, we 
would find that the ratio of rates 



k(b) 



k{b = 0) 



f Mb \ * 



Mb 2 



(21) 



where Q is a function of T/T G on the uncorrelated 
landscape, and on the correlated landscape is a func- 
tion of both T/T G and structural entropic factors hav- 
ing to do with the density of states of given similarity 
to a trap. So by inspection of eq. (pl|), in this low tem- 
perature regime the rate may go up or may go down 
with non-native interaction strength. 

However an important result arising from energetic 
correlations in the landscape is the existence of a crit- 
ical temperature T A where the dynamics becomes un- 
activated |39| . Above this temperature the dynamics 
is similar to reconfigurations in a normal liquid rather 
than the hopping dynamics of trap escape in a su- 
percooled liquid or glassy system (l^, |l8|, ^, [39| , i.e. 
above T A , the prefactor r _1 (6) remains nearly con- 
stant with ruggedness, since at these temperatures 
the Rouse modes depend much more weakly on the 
ruggedness introduced. The existence of such a tem- 
perature scale can be seen from the following simple 



argument []39| . We can think of escape from a trap as a 
mini- unfolding event: escape is driven by entropy and 
is opposed by the putative trap's low energy, say 2%. 
Then, as in unfolding, the escape barrier arises from 
a mismatch between entropy gains and energy losses 
as the system reconfigures out of the trap, so we can 
rewrite equation ( |lc| ) for the free energy relative to 
the state i as 

F(q) = qEi -TS c (q)-^(l- q) (22) 

where E N in (|lc| ) is replaced by £",;, Q in dig) is re- 
placed by q, defined as the fraction of contacts shared 
with state i, and density changes during untrapping 
are not particularly important since T A < Tg (ener- 
getic trapping occurs only when the polymer is col- 
lapsed |4j). The transition to unactivated dynamics 
occurs when the states typically occupied at that tem- 
perature have zero escape barrier. Setting E t in (22) 
to the thermal energy of states at temperature T,TT 



E, 



E(T) = - 



T 



(23) 



we note that T A occurs when the free energy profile 
is downhill away from the trap at q = 1, i.c when 
dFy/dQ = in our model, where the subscript > 
indicates the high q region of the piecewise free energy 
function (^2|) (eq. (||) for the configurational entropy 
has a piecewise structure). Using equations (m), (B2|), 
and (pJ3[), this gives 



T A = T n 1 



2fi) 



-1/2 



(24) 



for the transition temperature to activated dynamics. 
From ( p4|) we see that T A > T G by an amount which 
depends on the deviation from linear entropy loss over 
the total unconstrained entropy, i.e. by the entropic 
contribution to the barrier. There is no energetic con- 
tribution since we have used q as the order parameter 
and assumed pairwise interactions, p 



3 c.f. eq. ( |la| ) at Q = 0. For strata of states with Q > 
the larger ruggedness scale A(Q) increases T A and T G for 
that stratum of states. 

4 One must be consistent in interpreting eq. rt24|). In 
mean-field theory, <f> is extensive and T A > Tg in the ther- 
modynamic limit. But in the capillarity theory the en- 
tropic deviation cf> comes from surface entropy and should 
scale as N 2 ' 3 (or even with a smaller power if the inter- 
face is roughened). One might argue that since S ~ N, 
T A approaches T G as N — > oo, but matching the theories 
in this fashion is incorrect since eq. (E2I) is not valid in 
the capillarity limit. In the capillarity theory a dynami- 
cal transition can only be seen by investigating where the 
intensive surface tension vanishes. 



For typical size proteins of N ~ 100, T A w (1.6 — 
1.8) T G . A plot of the escape time on a correlated 
landscape is given below (see fig. ||). In the regime we 
are interested in, T G /T F « 1 (c.f. eq. (||)), so it is 
also true that 



!«' 



(25) 



and so the characteristic temperatures where folding 
occurs (~ T F ) are way above the transition temper- 
ature for activated diffusion by construction of the 
problem; see figure (p[). 

Expanding eq. (|Tj]) around 6 = using eq. ( |l0| ) 
gives the rate for weak non-native heterogeneity: 



k 



Sl + M Arf 



4z 2 



(26) 



where Ar/^ = rj{Q^) — ri{Q) and we've used Q^ = 1/2. 
The increase in rate occurs until around b A — 5b(N), 
where b A is where T F (b) = T A (b) (about 0.3 here), and 
6b(N) is the finite-size fluctuation of b A due to tem- 
perature fluctuations which round the transition |19| ] . 
That is, 



6b ^6T A = £- 



T 2 



Ta 



VM62 



Ta 



-7= (27) 



0. 



where eq. ( |lb| ) was used for the entropy at Q 
This gives a value for 8b w 0.02 - 0.04. 

Realistic values of b 2 /e 2 for a typical protein can be 
obtained from the ratio of folding to glass tempera- 
ture, given by 



T P /T G = A + y/X 2 -! 



(28) 



where A = y / z/2s (e/b) fJ5|| . Commonly accepted 
values of T F /T G for proteins are about 1.6 — 2.0. Using 
the values in table I for s and z gives 






proteins 



(0.1 - 0.15) 



(29) 



which is above the rate enhancement regime (see 
fig. pj) . If the effects of non- native rate enhancement 
are observable, they will be seen possibly in only the 
fastest folding proteins. On the other hand such an 
observation would support the existence of a dynamic 
glass transition in protein systems. 



IV. CONCLUSION 

As non-native heterogeneity is increased from zero, 
the folding rate initially increases by a factor of about 
2 — 4, then eventually drops rapidly towards zero (see 
fig. |l|). There is a regime near the Go Hamiltonian 



where the rate is relatively robust to changes in non- 
native interaction strength. The density of non-native 
polymer must be greater at the barrier peak than in 
the unfolded state for rate enhancement to be ob- 
served. Why does the rate initially increase? The 
upshot is as follows. First, it follows from energy 
landscape theory that if there is no change in density 
upon folding, then changes in ruggedness do not af- 
fect the barrier height for a well-designed protein (c.f. 
eq. (Jig)), that is, effects on the barrier must come from 
the coupling of non-native density to the degree of na- 
tiveness. Moreover, the strength of the ruggedness per 
residue increases with non-native density since there 
are simply more interactions, and non- native density 
tends to increase with nativeness, at least for weak 
non-native interaction strength. Then, since rugged- 
ness lowers the free energy, the free energy at the bar- 
rier position is lowered more than the unfolded free 
energy. So as non-native interaction strength is in- 
creased from zero, the barrier lowers and the rate in- 
creases if the effect on the prcfactor is weaker. But the 
prefactor is related to the reconfigurational diffusion 
time P,p5Lp39|, and since a dynamic glass transition 
is expected in such systems, there will be a window 
for weak ruggedness within which the diffusion time is 
relatively unaffected as ruggedness is increased. Thus 
the rate initially increases, as shown in fig. |l|. This 
phenomena provides a good example of how energy 
landscape theory can be applied to the physics of pro- 
tein folding to reveal and explain a counter-intuitive 
result. 
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V. APPENDIX A: CAVEATS DUE TO 

FINITE-SIZE, GENERIC ATTRACTION, AND 

STIFFNESS EFFECTS 

The derivation leading to eq. (M) assumed mean- 
field theory could be applied, that the protein could 
could be treated as a bulk system, and that proper- 
ties arising from chain connectivity would not alter the 
results arising from the energetics in the problem. In 
particular we have assumed that the polymer persis- 
tence length or Kuhn length £ K is much less than the 
length of a typical piece of disordered protein £ near 



the barrier peak, so that a dangling piece of disordered 
polymer may interact with itself. If the protein un- 
der study is particularly stiff and/or small, the return 
length of polymer fragments may be comparable to 
the length of the disordered pieces, reducing the num- 
ber of non-native interactions near the barrier peak 
relative to the number in the unfolded state. Then in 
eq. (||) the density rj{Q^) £ 7/(0) and no reduction in 
barrier height with non-native interactions would be 
seen. 

For the 27-mer lattice model £ K w 3 — 4; these mod- 
els are relatively stiff compared to their total length. 
For typical off-lattice models on the other hand, £ K « 
3 but they are considerably longer, e.g. for SH3, 
N = 57. If all the non-native polymer is in one strand, 
£ « N/2 at the barrier peak. Then £ K /£ « 0.26 
for the lattice model, and £ K /£ ~ 0.11 for the SH3 
off-lattice model. If the non-native polymer is dis- 
tributed among a number of disordered strands that 
can dress the native core, roughly (1/6) x N 2 / 3 |ll|, 
then £ K /£ « 0.40 for the lattice model, and £ K /£ ~ 
0.26 for the off-lattice model. 

Neither of these numbers are very small indicat- 
ing that the collapse transitions are quite rounded, 
and the effect on folding rate will be mild if it ex- 
ists. In fact the discrepancy of £ K /£ between the 
off-lattice and on-lattice models, although fairly small 
here, leads to different behavior of the non-native 
density rj(Q), as shown in fig. [|. In the lattice 
system rj{Q^) £ rj(Q v ), but in the off-lattice sys- 
tem rj(Q^) £ Tj{Qu)- Hence we anticipate the rate- 
enhancement effects will be seen in off-lattice mod- 
els, but probably not in at least the shorter on-lattice 
models B. 

Real proteins may tend to have some net homopoly- 
mer attraction inducing generic collapse. This de- 
creases the change in density upon folding and would 
further attenuate any rate enhancement effect present. 
However at least some proteins are sufficiently stable 
that collapse and folding are concomitant P7[| . More- 
over Go models, for which collapse and folding are 
concomitant by construction, give reasonably accu- 
ratepredictions of 0- values and barrier heights pl,[|,l7], 
|2[|21, 33, 34|. Collapse accompanies folding when the 
folding transition temperature T F given by eq. ( pn i s 
comparable to the collapse temperature Tg in eq. (QJ) . 
For weak ruggedness, 



To 



as n \ 2z e 2 2z a 2 



(30) 



So the effects of generic collapse are not important in 
the problem as long as e it as , to the first approxi- 
mation. 

Several additional features may affect the folding- 
rate. In finite-sized systems the unfolded state tends 
to have partial order. Moreover its position may drift, 
along with that of the transition state, as non-native 



variance is increased. This modifies the barrier height. 
Additionally the density r](Q) in figure [| is exact in 
the limit b — > 0, but for nonzero b, r](Q) may begin 
to alter in structure. Accounting for these effects will 
modify the folding rate, but shouldn't alter the general 
trend of figure H]. 



TABLE I. PARAMETERS IN THE MODEL 
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A. Figure Captions 



FIG. p] Rate vs. non-native heterogeneity is split up into two regimes, one where it assists folding, the 
other where it hinders. The rate plotted here is the folding rate at T F which is itself a function of b (c.f. 
eq. (pi), however T F changes by only ~ 1% over the range of this plot, (a) Schematic depicting the two regimes, 

(b) Result of the theoretical model introduced in the text (see equations [f(J, |2C], and figure |g) for a system 
with parameters given in table I. Initially the rate rises as ~ Nb 2 /e 2 , then strongly decreases for larger b as 
non-native interactions slow conformational transitions. Realistic protein interactions are believed to have 
typical values of b 2 /e 2 sa 0.1 — 0.15 [p5|, which is above the rate enhancement regime. The inset of (B) shows 
a semi-log plot of the same rate; it can be seen that there is a regime where the rate is roughly constant as 
ruggedness is increased from zero, then a turnover where the rate drastically decreases. 

FIG. H Energy, entropy, and free energy vs. Q in the model, for the Go model with A 2 = (solid line), and 
when ruggedness is introduced, when A 2 > (dashed line). We took b 2 = 0.04 (see eq. ||), where the folding 
rate is maximal (see fig. Ef) . Parameters used in the model are given in table 01 A bilinear approximation for 
the configurational entropy is used here, giving a tent functional form for the Go free energy at the transition 
temperature T F . The non-native density function r)(Q) used here is a fit to the off-lattice data in fig. 0C. The 
folding free energy barrier at T F is lowered by non-native heterogeneity, because the energy is lowered twice 
as much as the entropy (see eq. 1). 

FIG. [| As folding progresses, the non-native polymer halo surrounding the native core (central shaded 
globule) has more topological constraints placed upon it. Therefore the non-native packing density, given by 
the total number of non-native contacts divided by the characteristic volume of non-native polymer (open 
spheres), tends to increase. 

FIG. H Free energy profiles F(Q) and non-native polymer density r](Q) in the model for (a) the simple 
bulk-mean-field model used in the derivation, (b) a 27-mer lattice model, (c) an off-lattice model for the 57 
residue fragment corresponding to the a-spectrin SH3 domain (PDB code 1BK2). (a) Collapse occurs at Qe 
before the barrier peak at Q^. The transition is rounded for typical sized proteins, as in (b) and (c). In 
(b), the non-native density is overall fairly small, and is comparable in the unfolded and transition state, 
w(Qv) ~ viQ^) ~ 0.22. Thus eq. (Q) gives a barrier height roughly independent of b at least for small b. In 

(c) on the other hand, the non-native density rises to values larger in overall magnitude, and is monotonically 
increasing until the barrier peak: rj(Q^) w 0.35 and T](Q V ) ~ 0.2. For the parameter values in table I, eq. (|[) 
then gives a barrier height decreasing with b as AF^/T F m —10b 2 /e 2 . The drop-off at high Q in (B) and (C) 
is most probably due to stiffness effects on the small pieces of residual non-native polymer in this regime. 

FIG. |5j Log of the reconfiguration time vs. reciprocal temperature in units of T G , for a system of size 
N = 64, adapted from reference Q. This is used in equation ( p0| ) to produce the rate curve in figure |l|. 
At temperatures above T A the time to reconfigure is ~ r , below T A the time increases exponentially as 
exp(f(T/T A )N) with f(x) = for x < 1. The width of the transition — > as N — > oo and the value 
of T A — ►» 1.8 as JV — * co for the mean field correlated landscape. When T ;$ T G for a finite system, the 
deepest trap tends to dominate the kinetics, and the relaxation rate turns over to an Arrhenius law with slope 
corresponding to the barrier height for escape from that trap (shown schematically by the dashed line). 
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